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Abstract 

It is shown that monic orthogonal polynomials on the unit circle are the cha- 
racteristic polynomials of certain five-diagonal matrices depending on the Schur 
parameters. This result is achieved through the study of orthogonal Laurent poly- 
nomials on the unit circle. More precisely, it is a consequence of the five term 
recurrence relation obtained for these orthogonal Laurent polynomials, and the one 
to one correspondence established between them and the orthogonal polynomials 
on the unit circle. As an application, some results relating the behavior of the zeros 
of orthogonal polynomials and the location of Schur parameters are obtained. 
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I. Introduction 

It is very well known that orthogonal polynomials on the real line are the cha- 
racteristic polynomials of the principal submatrices of certain tri-diagonal infinite 
matrix called Jacobi matrix [1, 10, 19, 53]. The Jacobi matrix is just the represen- 
tation of the multiplication operator in the linear space of real polynomials when 
orthogonal polynomials are chosen as a basis. The fact that such polynomials satis- 
fy a three term recurrence relation is the reason for the tri-diagonality property of 
the Jacobi matrix. This property makes the Jacobi matrix simple enough to obtain 
properties of the orthogonal polynomials using operator theory techniques [2, 10, 
14, 20, 22, 29, 42, 43, 46, 53]. For example, one can get results about their zeros 
throughout the spectral analysis of the principal submatrices of the Jacobi matrix 
[6, 7, 11, 27, 28, 30, 31, 32, 33, 45, 52]. 

In the last years, there has been an increasing interest in the zeros of orthogonal 
polynomials on the unit circle due to their applications in discrete systems analysis, 
in particular, in digital signal processing [36, 38, 48, 49]. Unfortunately, only few 
things arc known about these zeros [3, 4, 5, 15, 17, 44, 47, 50, 51], because the situ- 
ation in the unit circle is rather more complicated than in the real line. Orthogonal 
polynomials on the unit circle satisfy too a three term recurrence relation that can 
be related to a tri-diagonal operator [5, 12, 13], but it does not provide a spectral 
representation for their zeros. However, it is possible to reach such a representation 
just by computing the matrix corresponding to the multiplication operator in the 
linear space of complex polynomials when orthogonal polynomials are chosen as a 
basis. The result is an irreducible Hessenberg matrix [19], much more complicated 
than the Jacobi matrix on the real line. So, this does not seem a so easy way to 
study properties of orthogonal polynomials on the unit circle. 

The aim of this paper is to improve this situation giving a five-diagonal matrix 
representation of orthogonal polynomials on the unit circle, which yields a spectral 
interpretation for their zeros. As we will see, this result comes from the matrix 
representation of the multiplication operator in the linear space of Laurent polyno- 
mials, when a suitable basis related to the orthogonal polynomials is chosen. This 
matrix representation gives a spectral interpretation for the zeros of orthogonal 
polynomials which is much simpler than the one given by the Hessenberg matrices. 
So, it provides an easier way to calculate these zeros and to study their behavior 
just using standard methods for eigenvalue problems of banded matrices. 

First of all we will fix some notations. For an arbitrary finite or infinite matrix 
M, M T is the transpose matrix of M, and M* — M . When M is a square matrix, 
M n means the principal submatrix of M of order n and, as usual, if M is finite, 
detM is its determinant. 

In what follows P :— C[z] is the complex vector space of polynomials in the vari- 
able z with complex coefficients. For n > 0, P„ :— (1,2,..., z n ) is the corresponding 
vector subspace of polynomials with degree less or equal than n, while P_i := {0} 
is the trivial subspace. As usual, if p € P„\P„_i, p* is its reversed polynomial, 
defined by p*(z) := z n p(z^ 1 ). A := C[z,z _1 ] denotes the complex vector space of 
Laurent polynomials (L-polynomials) in the variable z. For m < n, we define the 
vector subspace A TOi „ := (z m , z m+1 , . . . , z n ). Also, for any L-polynomial /, we will 
consider its substar conjugate defined by f*(z) = T := {z e C||z| = 1} and 

B> := {z e C||z| < 1} are called, respectively, the unit circle and the open unit disk 
on the complex plane. 
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Any hermitian linear functional £ on A (C[z n ] = C[z n ], n = 0, 1,2, . . . ) defines 
a sesquilinear functional (-,-)£:AxA^Cby 

(f,g)c := C[f(z)g(z-% /, g e A, 

and we say that /, g € A are orthogonal with respect to £ if (f,g)c = 0. The 
hermitian functional L is quasi-definite if there exists a sequence of orthogonal 
polynomials with respect to £, that is, a sequence (p n ) n >o in P satisfying 

(I) Pa € P n \Pn-l, 
(II) (p n ,Pm)c = £nS n , m , 4 7^ 0. 

The last condition can be replaced equivalently by 
(III) (p„,z fc ) £ = 0if0<fc<n-l, 
(Pn,* n )/:^0. 

Positive definite hermitian functionals (£ n > for all n) coincide with those 
given by C[f] := J c fdfj,, where /i is a positive measure with an infinite support 
lying on T. Due to this reason, the sequence (p n )n>o is called a sequence of orthog- 
onal polynomials on the unit circle [54] even in the general quasi-definite case. In 
particular, when £ n = {p n ,Pn)c — ±1 for all n, we say that [p n )n>o is a sequence 
of orthonormal polynomials on the unit circle. 

Given a quasi-definite hermitian functional C on A, (<p n ) n >o denotes the unique 
sequence of orthonormal polynomials with positive leading coefficients, whereas 
{4>n)n>Q is the unique sequence of monic orthogonal polynomials. They are related 
by ip n = K n (j) n , where K n := \(4> n ,4>n)c\~ 1 ^ 2 - Thus, if sg(-) is the sign function, 
((fin, (p n )c = sg((0„, 4> n )c)- I n what follows we consider C normalized by the con- 
dition £[1] = 1, so, <po(z) — (j>o{z) — 1 and Ko = 1. 

It is well known that the sequence (4> n )n>o is determined by the so called Schur 
parameters a n := <^ n (0) through the forward recurrence relation 

d>o(z) = 1, 

(1.1) 

<j) n (z) = z<t> n -i(z) + a n (j)* n _ 1 {z), n>\. 
From (1.1) we find that 

(1.2) — — = l-|a n |, n>l. 

(<Pn-l,<Pn-l)C 

and, therefore, e n := {<p n ,<p n )c/(<Pn-i,<Pn-i)c = sg(l - |a«| 2 ). Hence, apart from 
the first Schur parameter, that is always ao = 1, in the quasi-definite (positive 
definite) case it must be \a n \ ^ 1 (\a n \ < 1) for n > 1. Moreover, any sequence 
(a«)neN with this property yields, through the recurrence (1.1), a sequence of monic 
orthogonal polynomials on T, and the associated normalized functional C is unique. 
The parameters {afc}£ =1 C C\T allows to construct the finite set of polynomials 
{^fc}fc=o> t na t IS called the finite segment of orthogonal polynomials associated to 
{a,k}k=i- Given a finite segment {4>k}kZo °f orthogonal polynomials, any polyno- 
mial with the form <j> n (z) = Z(f) n -i(z) +t0*_ 1 (z),< e C\T, is called an extension of 
{(f>k}kZo- Extensions of a finite segment are the possible candidates to enlarge it. 
Relation (1.2) yields 

n 

(1.3) K n = Y[ Pk 1 , n>l, 

k=l 
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where p n := |1 — la^l 2 ! 1 / 2 for n > 1. With this notation we have the following 
forward and backward recurrence relations for the orthonormal polynomials and 
the reversed ones, 



(1.4) 
(1.5) 
(1.6) 
(1.7) 



Zlfn-^Z) = Pnfn(z) ~ a n ¥>n-l( z )> ^ > 1, 
¥>n-l(*0 = Pn*Pn( z ) - a n ZiPn-l(z), U>1, 

ip n (z) = a n (p* n {z) + p n ztp n -i(z), n > 1, 

ip* n (z) = a n ip n (z)+p n ip* n _ 1 (z), n>\, 



being p n := e n p n for n > 1. 

The finite segment of orthogonal polynomials associated to {ak} k =i let us define 
the n-th kernel 



(1.8) 



K n (z,y) ~ y^e k ip k (z)tp k (y), 



k=0 



where e n :— ((p n ,fn)c, that is, e = 1 and e n = YYk=i £ k for n > 1. Using the 
recurrence relation, the n — 1-th kernel can be written equivalently as 



(1.9) 



K n -i(z,y) 



¥n{z)Vn{y) ~ <fin(z)<Pn(y) 



zy 



if zy ^ 1, 



k e„z (y' n {z)ip n (y) - {iPn)'{z)Vn(y)) . if Z V = L 



In what follows it will play an important role the multiplication operator on A, 
defined by 

LT: A — > A 

/(*)-*/(*) 

Since it lets P invariant, using the orthonormal polynomials (tp n ) n >o as a basis of 
P, it is possible to obtain the matrix representation of the restriction of II to P, 
giving the following irreducible Hessenberg matrix [3, 15, 17, 18, 55] 



(1.10) 



where 



(1.11) 



H = 



■do,o d 0A 
di,o di,2 

$2,0 ^2,1 ^2,2 ^2,3 



- a j 



d-n, 



fl n+lIlL J + lft' if J = 0, 1,.. 

-a n a n+ i, if j = n, 

pn+i, ifj = n + l. 



and (a n ) n >o are the Schur parameters associated to the orthogonal polynomials 
considered. 

The characteristic polynomial of the principal submatrix TL n of 7i of order n is 
the corresponding n-th monic orthogonal polynomial [3, 15, 17, 55]. So, the spectral 
analysis of H n can give relations between the zeros of orthogonal polynomials and 
the Schur parameters, but, the fact that Tt is a Hessenberg matrix, together with 
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the complicated dependence of its elements with respect to the Schur parameters, 
makes difficult this task. Moreover, the numerical computations of zeros of high 
degree orthogonal polynomials, useful, for example, for applications in digital signal 
processing, becomes a non trivial problem due to the Hessenberg structure of Ji. 

So, a natural question that arises is how to find better spectral representations for 
orthogonal polynomials on the unit circle (we mean with this, the identification of 
any orthogonal polynomial as a characteristic polynomial of a matrix) . It would be 
desirable a situation so similar as possible to the one on the real line, in particular, 
we can think on finding banded spectral representations for orthogonal polynomials 
on the unit circle. Moreover, if we want to use this representations to connect the 
behavior of zeros of orthogonal polynomials and Schur parameters, we would need a 
simple dependence of the elements of the corresponding matrices with respect to the 
Schur parameters. We will find five-diagonal spectral representations for orthogonal 
polynomials on the unit circle satisfying this requirement. Apart from the obvious 
advantages for the study of orthogonal polynomials, this result implies a reduction 
of the eigenvalue problem for certain Hessenberg matrices to the eigenvalue problem 
of a five-diagonal matrix. 

The main idea to reach above result is to search inside the matrix representations 
of the full operator II on A. This leads to the study of basis of A related to 
orthogonal polynomials with respect a quasi-definite hermitian functional. As we 
will see, the most natural choice are those basis constituted by Laurent polynomials 
which are orthogonal with respect to the same functional. 

2. Orthogonal Laurent polynomials on T 

Orthogonal L-polynomials on the real line appeared in the early eighties in con- 
nection with the theory of continued fractions and strong moment problems [39, 
40]. Their study, not only suffered a rapid development (for a survey, see [37]), but 
it was extended to an ampler context, leading to a general theory of rational orthog- 
onal functions (see [8] and references therein). This theory cover, as a particular 
case, the orthogonal polynomials on T. However, the singularities of this particular 
case are lost in such a general theory, and, as we will see, these particularities are 
just the reason of their utility in the study of orthogonal polynomials on T. 

On the real line, orthogonal L-polynomials are a natural generalization of the 
orthogonal polynomials when the related functional, initially defined only for poly- 
nomials, is extended to the space of L-polynomials. Our interest is in the general- 
ization of this idea to the unit circle, where the corresponding functional is already 
defined for the full space of L-polynomials. 

Although we will deal with orthogonal L-polynomials with respect to a general 
quasi-definite hermitian functional on T (see [23] for the analogous generalization 
on the real line), just to understand the following definition, let us consider a 
positive definite hermitian functional C on A. Then, the sesquilinear functional 
(■,■)•£ is an inner product on A, and the orthogonal polynomials with respect to C 
appear from the standard orthogonalization of the set {1, z, z 2 , . . . }. Analogously, 
using the Gram-Schmidt procedure we can get an orthogonal basis {/ ra }„>o of 
A starting from the ordered basis {1, z, z^ 1 , z 2 , z~ 2 , . . . } of A. If we define the 
subspaces A^ := A_„ ; „,A^ n+1 := A_ ni „ + i for n > and Aj^ := {0}, then such 
an orthogonal basis satisfies 
(la) /„ G A+\A+_ 15 
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(Ila) (/„, f m ) C = i n S n , m , In ^ 0. 

So natural than this is to start with the ordered basis {1, z~ x , z, z~ 2 , z 2 , . . . }. 
In this situation, the Gram-Schmidt orthogonalization process gives an orthogonal 
basis {f n }n>o 01 A satisfying 
(lb) /„ G A-\A"_ l5 

(lib) (/„, f m ) C = t n Sn, m , in + 0, 

where Aj„ := A_„_„, A 2rl+1 := A_„_i )n for n > and AZ X := {0}. 
Above discussion is the origin of the following definition. 

Definition 2.1. A sequence (f n ) n >o in A is called a sequence of right (left) or- 
thogonal L-polynomials on T if 

(I) fn G A+ ( - ) \A ? |i 1 ) 

and there exists a hermitian functional £ on A such that 

(H) {fn: fm)c — ^-n^n,m > ^n 7^ ^' 

Then, we say that (f n )n>o is a sequence of orthogonal L-polynomials with respect 
to C. If, in addition, l n = ±1 for all n, then we say that (f n )n>0 is a sequence of 
orthonormal L-polynomials. 

Remark 2.1. Similarly to orthogonal polynomials, Condition (II) in Definition 2.1 
can be replaced equivalcntly by 

(Ilia) (f 2n , z k ) c = if -n + 1 < k < n, (f 2n , z~ n ) c £ 0, 

(f2n+i,z k ) c = if -n < k < n, (/ 2 „+i,z n+1 ) £ ? 0, 
in the case of right orthogonal L-polynomials. For left orthogonal L-polynomials 
the equivalent condition is 
(IIIb) (/ 2n , z k ) c = if -n < k < n - 1, (f 2n , z n ) c ? 0, 

(f2n+i,z k )c = if -n < k < n, (f2n+i,z- n - 1 ) c 0. 
Analogously to orthogonal polynomials, right and left orthogonal (orthonormal) 
L-polynomials are unique up to non null (unimodular) factors. 

Contrary to what happens in the real line, in the unit circle right and left or- 
thogonal L-polynomials are closely related. 

Proposition 2.1. Let £ be a hermitian functional on A and let (f n ) n >0 be a 
sequence in A. Then, (f n )n>0 is a sequence of right orthogonal (orthonormal) 
L-polynomials with respect to C iff (f n *)n>o is a sequence of left orthogonal (or- 
thonormal) polynomials with respect to C. 

Proof. Obviously /„* G A^\A^_ 1 iff /„ G A+\A+_ 1 . The rest of the proof is just a 
consequence of the hermiticity of £, since it implies (/„*, f m *)c = (fn, fm)c- ^ 

Moreover, in the unit circle, orthogonal L-polynomials can be easily constructed 
from orthogonal polynomials, something that does not hold in the real line. This 
fact, although trivializes such rational functions, is the key for their usefulness in 
the study of orthogonal polynomials. 

Proposition 2.2. Let C be a hermitian functional on A and let (f n ) n >o be a 
sequence in A. Let us define 

PtnW - ^J^z- 1 ), pt n+1 (z) = z n f 2n+1 (z), n > 0, 

P2n( Z ) = z ' n f2n{z), P 2 n+l( Z ) = ^1 2n+\ ( z ~ 1 ) . " > 0- 
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Then, (f n )n>o is a sequence of right (left) orthogonal L-polynomials with respect to 
£■ iff {Pn^ ')n>o is a sequence of orthogonal polynomials with respect to C More- 
over, (f n )n>o are orthonormal iff (pi^ '*)n>o so are. 

Proof. It is straightforward to prove that p+ £ P„\P„_i iff /„ £ A+\A^_ X . More- 
over, (pj n ,z k ) c = (f 2n ,z n - k ) c and {pt n+l ,z k ) c = (f 2n +i, z k - n ) c . Thus, p+ n is 
orthogonal to 1, . . . , z 2 ™ -1 iff fm is orthogonal to z~ n+1 , . . . , z n , and pt n +i is or_ 
thogonal to 1, ... , z 2n iff f2 n +i is orthogonal to z~ n , . . . , z n . Besides, (p^,, z 2n )c = 
(fin,z- n ) c and {ptn+n z2n+1 ) c = (hn+i, z n+1 )c- Therefore, according to Remark 
2.1, (Pn )n>o is a sequence of orthogonal polynomials iff (f n )n>o is a sequence of 
right orthogonal L-polynomials. Finally, since (p+,p+)£ = (f n ,fn)ci we have that 
(Pn)n>o arc orthonormal iff (f n ) n >a so are. A similar proof works in the case of 
left orthogonal L-polynomials. □ 

Remark 2. 2. Above result establishes in the unit circle a one to one correspondence 
between sequences of orthogonal (orthonormal) polynomials and sequences of right 
or left orthogonal (orthonormal) L-polynomials. So, any sequence of right orthogo- 
nal (orthonormal) L-polynomials (f n ) n >o is obtained from a sequence of orthogonal 
(orthonormal) polynomials (p n )n>o by the relations 

f 2n (z) = z- n p* 2n (z), n>0, 

/2n+l(z) = Z~ n p 2n +l(z), Th > 0. 

The corresponding sequence of left orthogonal (orthonormal) polynomials (f n *)n>o 
is given by 

f2n*(z) = Z~ n p 2n (z), n > 0, 

hn+i*{z) = z- n - 1 p* 2n+1 (z), n>0. 

Moreover, £„ := (/„, f n ) c = (/„*, f n *)c = {Pn,Pn)c- 

Notice that /„ and /„* are linearly independent for n > 1: for odd index it is 
obvious; for even index it is a consequence of the same property for p n and p* . 

From above comments we see that the conditions for the existence of orthogonal 
polynomials or L-polynomials are the same, that is, quasi-dehnite hermitian func- 
tionals on A are just those hermitian functional for which there exist (right or left) 
orthogonal L-polynomials. 

Definition 2.2. Given a quasi- definite hermitian functional £ on A, we denote by 
(Xn)n>o the sequence of right orthonormal L-polynomial defined by 

X2n(z) := z~ n ip* 2n (z), n>0, 

X2n+l{z) := Z~ n ip2n+l{z), U > 0, 

where (f n )n>o is the corresponding sequence of orthonormal polynomials with posi- 
tive leading coefficients. We refer to (xn)n>o ((Xn*)n>o) as the standard right (left) 
orthogonal L-polynomials associated to C 

The standard right and left orthonormal L-polynomials satisfy some useful rela- 
tions that are direct consequences of the recurrence relation for the corresponding 
orthonormal polynomials. 
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Proposition 2.3. Let C be a quasi-definite hermitian functional on A and let 
(Xn)n>o be the related sequence of standard right orthonormal L-polynomials. Then, 

(i) ( X2n{ ?\ ) = — h n 1 ) ( X2n ~ l{ ?\ ), n > 1, 

\X2n*{z) J p 2n \ 1 «2n / \ X2n-1* (z) J 

( X2n -}^)=e 2n ( X ^-^), n>l, 

V X2n(z) J \ X2n*(z) J 

(Hi) z( ) = e 2n+1 ( X ^?A, n>0, 

\X2n+l*{z) ) ^ \Xln+\(z) ) 

6 n := ( M , n> 1, 

where (a n ) n eN are the Schur parameters associated to L, p n = |1 — lanpl 1 / 2 and 

Pn = £npn with £ n = Sg(l - |ft„| 2 ) . 

Proof. All the relations follow straightforward from (1.4), (1.5), (1.7) and Definition 
2.2. ' □ 

2.1. Recurrence relation for orthogonal L-polynomials on T. 

We know that, in the unit circle, the action of the multiplication operator over the 
orthogonal polynomials does not provide a recurrence relation for them. However, 
as we see in the next proposition, for the orthogonal L-polynomials the situation is 
much better. 

Proposition 2.4. Let (f n ) n >o be a sequence of right orthogonal L-polynomials on 
T. Then, there exist n n ^ k GC, \n — k\ < 2, such that 

n+2 

Zfn(z)= ^ Tfn,kfk(z), n>0, 
k=n-2 
n+2 g 

zf n *(z)= nk, n fk*(z), n>0, 

k=n— 2 

where £„ — (/„, f n )c and we use the convention f k = for k < 0. 

Proof. Notice that zf n G zA+ C A+ +2 = (/ , /i, . . . , /„+2>. Moreover, since /„ is 
orthogonal to (/o, /i, . . . , fn-i) = A n _i, we find that zf n is orthogonal to zA rl _i D 
A„_3 = (fa, fi, ■ ■ ■ , /n-3>- Therefore, zf n e (f n -2, ■ /n+2). Taking into account 
that (f n *)n>o is a sequence of left orthogonal polynomials, we find following similar 
arguments that zf m e (fn-2*, ■ ■ ■ , fn+2*)- Therefore, for n > 0, 

n+2 n+2 

Z.fn(z) = ^ 7r n,kfk(z), Zf„*(z) = ^ ^n,k*fk*(z), 
k—n — 2 k—n — 2 

with 7r„, fe = (zfn, fk)c/(fk, fk)c and n n ,k* = (zf n *,fk*)c/(fk*,fk*)c- Using the 
definition of the substar conjugate, we find that 7r„ ; fc* — ftk,n(-n/ (-k- □ 

Above proposition says that orthogonal L-polynomials on the unit circle satisfy 
a five term recurrence relation. This fact was already known for orthogonal L- 
polynomials on the real line, and used to solve the strong Hamburger moment 
problem through operator theory techniques [21]. 
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From the relation between orthogonal L-polynomials and orthogonal polynomials 
on T, it is possible to obtain explicitly the coefficients of the recurrence relation for 
the orthogonal L-polynomials in terms of the Schur parameters associated to the 
orthogonal polynomials. 

Proposition 2.5. Given a quasi- definite hermitian functional L on A, the related 
standard right orthogonal L-polynomials (Xn)n>o satisfy the recurrence relation 

zxo(z) = -aiXoO) + PiXi(z), 

z( x2n -^)=m n -J X2n -f\)+M 2n ( y n >i, 

V X2n{z) J \X2n-l(z) J \X2n+l{z)J 
( —P n a n +1 PnPn+l \ M ( ~Pn a n+l PnPn+l 

™ \ —a n a n+ i a n p n +i J ' ™ \ —a n a n+ i a n p n+ i 

where (a„)„ eN are the Schur parameters associated to L, p n = |1 — |a Tl | 2 1 and 
Pn = £nPn with s n = sg(l - |a„| 2 ). 

Proof. From equations (1.4), (1.7) and Definition 2.2, we have that, for n > 1, 

ZX2n-l = Z 2 ~ n <p 2 n-1 = Z 1 ^ n (p 2 nV2n - a 2 nV 2 n-l) 

= Z n P2n(P2n+l ( P2n+l — a 2ri+l</?2ri) — Z 1 ™a 2 „ (a 2n - 1 l^2n-l + p2ri-l</?2n-2) 
= p2nP2n+lX2n+l — P2nO<2n+\X2n ~ ~Q>2n-\&2nX2n-\ — p2n-lO-2nX2n-2, 
ZX2n = Z X ~ n <p 2n = Z 1 ~ n {a 2n ^2n + P2n^>2n-l) 

= Z 0,2n\P2n+W2n+l - a2n+W2n) + z P2n{a2n-l<P2n-l + P2n-l¥>2n-2) 
= a 2n P2n+lX2n+l — a 2n a 2 n+lX2n + a 2n -ip 2n X2n-l + p2n-lf>2nX2n-2- 

Besides, from Definition 2.2 and (1.4), 

zxo = z<fo = zipo = pupx - a x ipl = pxx\ ~ aiXo, 
which completes the proof. □ 

Remark 2.3. Following similar arguments to previous proof, from (1.5) and (1.6) 
we find that the recurrence relation for the standard left orthogonal L-polynomials 

(Xn*)«>0 IS 

\Xi*(z)J \ Pi J \X2*{z) J 

z( ^ ( ^)=MjJ x2 "-y^v^« + if X2n+i 1 z !y »>i- 

\X2n+l*{z) J Zn \ X2n*(Z) J \X2n+2*{z)J 

3. Orthogonal polynomials on T and five-diagonal matrices 

The recurrence relation for the standard orthonormal L-polynomials provides a 
five-diagonal infinite matrix that plays in the unit circle a similar role to the one 
played by the Jacobi matrix in the real line. 



n > 1, 
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Definition 3.1. The five-diagonal matrix T associated to a quasi- definite hermit- 
ian functional £ on A is the following infinite matrix 





M 2 



( -a\ Pi 
M? 


V : : 
/ -oi Pi 
-pia 2 — aia 2 
pifa aip 2 


V 



...X 






M 4 


-P2C13 P2P3 
-a 2 a 3 a 2 p 3 







-a3 ft 4 -P4«5 P4P5 

fa Pi a 3 fa -a 4 a 5 a 4 p 5 



where (o„)„ eN are t/ie Schur parameters related to C, p n = |1— |a„| 2 | and p n = e n p n 
with e n = sg(l - |a„| 2 ). 

Remark 3.1. Proposition 2.5 means that T is just the matrix of il with respect 
to the basis of A constituted by the standard right orthonormal L-polynomials 
(Xn)n>o related to C. When C is positive definite, Proposition 2.4 implies that T 7 " 
is the matrix of II when the standard left orthonormal L-polynomials (xn*) n >o are 
the basis chosen for A. Taking into account this proposition we get that, in the 
general quasi-definite case, the matrix of II with respect to (xn*)«>o is given by 
= ET T E 1 being E the infinite diagonal matrix 



E := 



/e 

ei 

e 2 

V ; ; ; 



with e„ = (x«, Xn)c — {<Pn> <Pn)c- In fact, from Remark 2.3, we find that T* is the 
result of substituting in T T the coefficients p n by p n and vice versa. 

Notice that II = IlWn( 2 ) = LT^nW, where IlW, i = 1,2,3, are the linear 
operators on A defined by 



n^: A — >A , 



II< 2 >: A^A , 



II< 3 >: A^A , 



From this fact and Proposition 2.3, wc find that T = T^T 1 ^ and ;F* = T^T^ 
where T^ x \ are the following block-diagonal matrices 



:= 



/0i 
e 3 ••• 
e 5 ■■• 



:= 



(\ 

e 2 

e 4 

V: 



■■J 



being all the blocks two-dimensional, excepting the first one for fW, that is one- 
dimcnsional. 
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This gives a decomposition of the five-diagonal matrices T, T* as a product of 
two tri-diagonal matrices with special properties: since 9 n 0„ = I2 for all n, we 
have that T^T^ = I = ~T^T {i) , i = 1,2 (7 is the infinite unit matrix). In 
the positive definite case p n = p n and, thus, T^\ T^ are symmetric and, so, 
unitary too. Therefore, T is unitary for a positive definite functional C. This is 
not a casuality, since, if p is the measure on T associated to the positive definite 
functional £, then T is the matrix of the unitary operator 

with respect to the Hilbert basis (xn)n>o of the space of /z-square integrable func- 
tions L 2 n . 

Notice that, for the principal submatrices of order n, we can write too T n = 

T^Tn^ and T m — Tn^Tn\ but now wc can only state that Tn" 1 Tn"' — I n for 

even n and Tn^Tn^ = I n for odd n. 

Analogously to what happens for orthogonal polynomials on T and the cor- 
responding Hessenberg matrix, one would expect a relation between the zeros 
of orthogonal L-polynomials on T and the principal submatrices of the related 
five-diagonal matrix. But, the connection between orthogonal polynomials and L- 
polynomials on T, provides finally a relation of those principal submatrices with the 
zeros of orthogonal polynomials. This fact justifies the mentioned analogy between 
the Jacobi matrix on the real line and the the five-diagonal matrix found on the 
unit circle. 

Theorem 3.1. Let T be the five-diagonal matrix associated to a quasi-definite 
hermitian functional on A with monic orthogonal polynomials (4> n )n>a o,nd standard 
right orthonormal L-polynomials (Xn)n>o- Then, for n > 1: 
(i) The characteristic polynomial of the principal submatrix T n of T of order n 
is (j) n . 

(ii) The eigenvalues of T n have always geometric multiplicity equal to 1. 

(iii) An eigenvector of T n corresponding to the eigenvalue A is given by A„(A), 
where X n (z) := z^\xo(z),Xi(z), ■ ■ ■iXn-i{z)) T . 

Proof. From Proposition 2.5, and using 2.3 (iii), we can write 



(zL n - T n )X n (z) = b n (z), 

PnZ 1+[s ~^ ] Xn*{ z ) (0)°) - - - , 0, 1) T , if n is even, 
PnZ 2 Xn (z) (0,0,..., 

Pn — 1 1 &n—l 

) T , ifnisodd, 



b n (z) 



where /„ is the unit matrix of order n. 

If n is even, applying Cramer's rule to solve above system with respect to 
Xn-i(z), we get 



Xn-i(z) = - / — ^det 
det(z/„ - T n ) 



( \ 

zl n -l — T n -1 



V ••• PnZXn*(z)J 



det(z/»_i - T n -i) . . 

1 J i-l t t- \ PnZXn*\ z ). 

det(z/„ - Tn) 
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Since ifj(z) = Kj(f)j(z) and pj — kj-i/kj, using Definition 2.2 we find that 

4>n{z) = 4>n-l{z) 

det(2;/„ - T n ) det(z/„_i - T n -\) ' 

When n is odd, Cramer's rule to solve the initial system with respect to \n-i{z) 
gives 

/ \ 

zI n -2 - T n -2 ■ ■ 



Pn-lPnXn(z) Pn-\a n 

V ••• a n -ip n Xn{z) z + a n -ia n J 



det(z/„ - T n ) 



det(zI n - 2 - F n -i) , , 

ZPn-lPnXn{z), 



det(z/„ - T n ) 
and, therefore, 



(fin(z) 4>n-2(z) 



det(zl n - T n ) det(z/„_ 2 - T n - 2 ) 
Hence, we find by induction that, for n > 1, 

4>n(z) (f>l{z) 



det(z/„ - T n ) det(z/i - T\) ' 

which proves (i) since <fii(z) = z + a\ = dct(z/ 1 — Ti). 

So, the eigenvalues of T n coincide with the zeros of </>„. Using Definition 2.2, we 
find that 



bn(z) 



p n ip n (z) (0, 0, ... , 0, 1) T , if n is even, 

p n (p n (z) (0,0, . . . ,p„-i,a„-i) T , if n is odd. 



Thus, if A is an eigenvalue of T n , it must be (XI n — T n )X n (X) — and, hence, 
just showing that A"„(A) ^ 0, (iii) is proved. First of all, notice that X n (z) is well 
defined for any value of z since their components are polynomials in z. If A ^ 0, 
the first component of X n (X) is non null because Xo( z ) = 1. On the contrary, when 
A = 0, one of the last two components can not vanish because 



X n (0) = 



(0, 0, . . . , K n -2, K„-ia„_i) T , if n is even, 
(0, 0, . . . , 0, k„_i) t , if n is odd. 



It only remains to prove (ii). If v = (vi, v 2 , ■ ■ ■ , v n ) T satisfies T n v = Xv, then 
-axvi + piv 2 = Xvi, 



( V2k - 1 ) + M 2k ( V ^)=x( ™ ), *=!,. ..,[ S ]-1, 

V V 2 k J \V 2 k+2/ \V 2k +l J 



-Pn-ia n v n -i - an-ia n v n = Xv n , if n is even, 

— T / v n - 2 \ f-p n -ia n \ .( v n -i\ . 

Mi_ 2 + I )v n = A I , if n is odd. 

\Vn-iJ \-a n -ia n J V v„ 



M.J. CANTERO, L. MORAL, L. VELAZQUEZ 



13 



Let us suppose first that A ^ 0. In this case we are going to prove by induction 
that i>i = iff v n = for all n. When v\ — 0, the first equation of previous system 
gives v 2 — 0, while the second equation implies that, if v 2k -i = v 2k = for some 
k < [n/2] - 1, then 

M 2k ( V2k+1 ) = x( ° V 

Using the expression of M n given in Proposition 2.5, we get from this identity that 
v 2 k+i = v 2 k+2 = 0. Therefore, if v\ — 0, then v 2 k-i = v 2k = for k < [n/2], what 
proves that v — if n is even. When n is odd, the last identity of the system for v 
gives v 2n+ i — and, so, v — too. 

Hence, if v = (vi, v 2 , . . . , v n ) T , v' = (v[, v' 2 , . . . , v' n ) T are eigenvectors of T n with 
the same eigenvalue A ^ 0, then v\,v[ ^ 0. Since w = v[v — viv' satisfies F n w = Xw 
and wi = 0, it must be w = and, thus, v, v' are linearly dependent. 

Let us consider now the case A = 0. If v = (vi, v 2 , ■ ■ ■ , v n ) T satisfies T n v = 0, 
then, taking into account that a n = 0, we have that 

' -aivi + p!v 2 = 0, 
Mi k _ 1 ( V2k - 1 ) + M 2 J V ^)=0, fc = l,...,[|]-l, 

< \ V 2k J V V 2 k+2 J 

Mj_ 2 ( Vn ~ 2 ) = 0, if n is odd. 

I \V„-lJ 



The last equation is equivalent to p n ~2 v n-2 + «n-2 !) n-i = 0, while, multiply- 
ing on the left by £>2k, the second relation becomes f>2k-iV2k-i + &2fc-i^2fc = 
— a 2 k+\v 2k+ i + p 2 k+iv 2 k+2 = 0. Therefore, the system for v is equivalent to 



e 2 k- 1 ( V2k - 1 )=o, fc = i,...,[^i], 
V v 2k j 

-a n -iv n -i + Pn-iVn = 0, if n is even. 



from what it is straightforward to see that v must be proportional to X n (0). □ 

Remark 3.2. From Remark 2.3, it can be shown that (i) and (ii) remains true for 
the matrix T*. Besides, the eigenvectors associated to an eigenvalue A of T m are 
given by X nt (X), being X w (z) := {xo*(z), Xi*(z), ■ ■ ■ , Xn-i*(z)) T ■ Notice that, 
when A = 0, 

x ^ _ | (0,0, . . . ,0,k„_i) t , if n is even, 

I (0, 0, . . . , K n - 2 , K. n -\a n -\) T , if n is odd. 



Remark 2.1 implies that T m = E n T^E n , where E n is the principal submatrix of 
E of order n. Thus, if A is an eigenvalue of T n , and taking into account that E\ = 1, 
we find that J^E n X n *(X) = XE n X n *(X), that is, E n X n *(X) is an eigenvector of 
T% with eigenvalue A. 

Notice that, if A is a zero of <j) n , then we can take as associated eigenvectors of 
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T n and F tn , V n (X) and V^*(A) respectively, where 



' (xo(z),Xi(z), ■ ■ ■ ,Xn-i(z)) T , ifA^O, 



V n (z) = < (0,0, . . . ,/9 n _i,a„_i) T , if A = 0, even n, 

k (0,0, ...,0,1) T , if A = 0, odd n, 

' (xo*(z),Xi*(z), ■ ■ ■ ,Xn-i.(z)) r , if A t^O, 



V n *(z)=l (0,0, ...,0,1) T , if A = 0, even n 

(0, 0, ... , p n -i, a„_i) T , if A = 0, odd n. 



Taking into account Remark 3.1, the eigenvalue problem for T n can be translated 
into 



That is, the zeros of <j> n can be viewed as the eigenvalues of the five-diagonal 
matrices T n and T mi or, alternatively, as the generalized eigenvalues of the tri- 



Previous theorem gives a spectral interpretation for the zeros of orthogonal poly- 
nomials on T, that allows to calculate them using eigenvalue techniques for banded 
matrices. This implies a reduction of their computational cost if compared with 
the calculation using Hessenberg matrices [19]. The banded structure of the five- 
diagonal matrices, together with their simple dependence on the Schur parameters, 
permits even to obtain properties of the zeros by means of standard matricial tech- 
niques, as we will show afterwards. In fact, this banded structure makes possible 
to apply similar techniques to those usual for the Jacobi tri-diagonal matrix on the 
real line. 

Besides, contrary to the Hessenberg matrix TL associated to the orthogonal poly- 
nomials on T, every Schur parameter appears in only finitely many elements of the 
matrix T . This makes easier for T than for H the analysis of the effects of pertur- 
bations of the sequence of Schur parameters. In particular, every modification of a 
finite number of Schur parameters induces a finite dimensional perturbation of the 
five-diagonal matrix T, something that is not true for the Hessenberg matrix TL. 
We will take advantage of these facts in the following section. 

As for the relation between T n and TL n , the geometric multiplicity of any eigen- 
value of an irreducible Hessenberg matrix is always one [19] and, so, it follows from 
Theorem 3.1 that T n and TL ni not only have the same characteristic polynomial, 
but are indeed similar matrices. This fact was not obvious since T n and TL n are 
matrix representations of different truncations of the multiplication operator. 



As a first application of the five-diagonal representation for orthogonal polyno- 
mials on T, we will derive bounds for their zeros in the general quasi-definite case, 
where only very few things are known. 

Theorem 4.1. Let L be a quasi- definite hermitian linear functional on A, (a n ) ne iq 
the corresponding sequence of Schur parameters and {z™}™ =1 the zeros of a n-th 
orthogonal polynomial associated with C. Then, 




diagonal pencil (fn^jfn^) or (Fn\ J 7 ^) depending on if n is odd or even. 



4. Applications 



Ri < h i < i?2, 1 < j < n => K x < \z] 



<K 2 , l<j<n, 
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where K x = Rf + R% - K 2 , K 2 = (R 2 + K) 2 , K = max { 1 1 — Rj | 1/2 , |1 - Rl\ 1/2 } ■ 

Proof. Applying Gcrshgorin theorem [19, 24] to the matrix T n we find that their 
eigenvalues have to lie on a union of disks Dj, j = 1,2, ... ,n, with centers 



f -oi if J = 1, 

\ -aj-xdj, if j = 2, . . 



and radii bounded by 



maxp fe , if i = 1, 

(max/} fc ) 2 + 2(maxp/ c )(max|a/ £ |), if j = 2, . . . ,n. 

k<n k<n k < n 

If i?i < |aj| < R 2 for 1 < j < n, then i?i < |ci| < R 2 , n < K, and R\ < \ Cj \ < R\, 
rj < K 2 + 2KR 2 , for 2 < j < n. Since \cj\ — rj < \z\ < \cj\ + rj for z e Dj, any 
eigenvalue A of JF„ must satisfy minji^!, Ri — K} < |A| < max{K 2 , R 2 + K}. The 
theorem follows from the fact that Ki < R\ — K and K 2 > R 2 + K: K 2 > R 2 + K 
\&R 2 +K>1, which is true since R 2 + K > R 2 + |1 - R 2 \ x l 2 >\;\iR 2 < \ then 
R x < 1 too and, thus, K x = R\ + R 2 - K 2 < R^ + R 2 - (R 2 +K) = R 1 -K; when 
R 2 > 1 we have that K x < Ri - K iff (2R 2 - 1)K > R{ - R 1 - K 2 , which is true 
since Rl~ R±- K 2 < Rj- Ri-ll- Rj\ < 0. □ 

So, bounds for the complete sequence of Schur parameters give uniform bounds 
for the zeros of orthogonal polynomials. Notice that, when applying Gershgorin 
theorem to the principal submatrices H n of the Hessenberg matrix H, we do not 
get in general uniform bounds for the zeros because each new row includes more 
non vanishing elements. In fact, using (1.10) and (1-11) we would have found that 
the centers of the corresponding Gcrshgorin disks are the same ones given before, 
but the bounds for the radii would be now 



( (maxp fc ) if j = 1, 

k<n 

(maxpfe) + V"(max / 9 fc ) r (max |a fc |) 2 , if j = 2, . . . ,n, 

k<n z — ' k<n k<n 

r=l 



which does not always give uniform bounds for the zeros in the quasi-definite case. 

The bounds that Theorem 4.1 gives for the zeros of orthogonal polynomials 
locate them in an annulus of radius K 2 — K\ = R 2 — R\ + 4KR 2 + 2K 2 . Thus, the 
best bounds appear when K is close to 0, that is, when R\,R 2 are close to 1. Since 
K\, K 2 — > 1 when R\,R 2 — ► 1, we have the following immediate consequence. 

Theorem 4.2. Let C be a quasi- definite hermitian linear functional on A, (a n ) ne ^ 
the corresponding sequence of Schur parameters and the zeros of a n-th 

orthogonal polynomial associated with L. Then, Ve > 0, 36 > such that 

\\aj\-l\<6, 1 < j < n ||z"| - l| < e, 1 < j < n. 



One choice that ensures this fact is 6 = y 1 + 4(i+ e ) — ^ • 

Proof. Following the notations given in Theorem 4.1, now Ri = 1 — 6, R 2 = 1 + 6, 
with 6 > 0, and, thus, K = ^J26 + 6 2 and K 2 = 1 + e, K x = 1 - e + 6 2 , where 
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e = 2K(1 + S + K). As K — > when 5 — ► 0, the implication in the theorem turns 
out to be true. The value given there for 5 is just the only positive solution for the 
equation 2^/25 + 6 2 (l + S + V28 + 5 2 ) = e. □ 

Roughly speaking, this last result says that, when the Schur parameters are close 
to the unit circle, the zeros of orthogonal polynomials so are. This fact is easy to 
derive for positive definite functional because, in this case, both, Schur parameters 
and zeros, lie on the open unit disk and, hence, the relation 



(-irn 



implies that 



< z'- 



< 1 for 1 < j < n. Theorem 4.2 is the generalization of 



this property for the quasi-definite case. 

In this first application we have exploited the banded structure of the five- 
diagonal matrix representation for orthogonal polynomials. Now we will show the 
advantages of this representation for the analysis of perturbations. 

Let us consider the five-diagonal matrix T of Definition 3.1 as a function of 
the sequence (a n )neN- Then, its principal matrix of order n becomes a function 
J>i(oi)tt2) • ■ • ,dn) of n complex variables. Given {aj}™ =1 C C\T, the analysis of 
the spectrum of T n {a\, 02, . . . , a n ) is equivalent to the study of the zeros of the last 
polynomial in the associated finite segment of orthogonal polynomials {</>j}™ =0 . 
We are going to analyze the variation of these zeros under perturbations of the 
parameters {oj}™ =1 . 

Let {aj(t)}% =1 C C\T be a set of parameters depending on a real or complex 
variable t in a neighborhood of t = 0. We consider aj = aj(0) as unperturbed 
parameters and use for them previous notations. For the perturbed ones {aj(£)}£ =1 
we adopt following notations: {(/>*• }™_ is the perturbed finite segment of orthogonal 
polynomials, the corresponding orthonormal polynomials are given by </jj (z) = 
Kj(t)<j>j(z) with Kj-i(t)/Kj(i) = pj(t) = |1 - Idj^)! 2 ! 1 / 2 , and K l n is the related 
n-th kernel. The associated standard right and left orthogonal L-polynomials arc 
denoted by \j an d Xj* respectively. Also, E n (t) is the diagonal matrix given by 



E n {t) 



( l 

ei(t) 



\0 



\ 



e n -i{t)/ 



fe=i 



where ej(t) = sg(l— \a,j(t)\ 2 ). Besides, we write F n (t) := T n (ai{t), a 2 {t) 1 . . . , a n {t)). 
Finally, and are the related vector polynomials given in Theorem 3.1 and 
Remark 3.2, which yields the eigenvectors of F n (t) and F m (t) := E n (t)F n (t) T E n (t) , 
respectively. The same holds for the vector polynomials V* and . 

Our aim is to study the variation of the zeros of 4>n as a function of the (real or 
complex) variable t for suitable perturbations. For this purpose it would be useful 
a "Hellmann-Feynman" type theorem [25, 26, 34, 35] for the eigenvalues of F n {t). 
Although F n (t) is not self-adjoint and, even, neither normal too, we can find such a 
result taking advantage of the relation between the eigenvectors of F n (t) and F n (t) T 
[19] given in Remark 3.2. This is the basic idea under the proof of the following 
proposition. 
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Proposition 4.1. If cij(t) and cij(t) are differentiable at t = for j = 1, 2, . . . , n, 

and we can express an eigenvalue of F n (t) := T n {ai{t),a2{t) 1 . . . , a n (t)) as a func- 
tion \{t) differentiable at t = 0, then 



V W {\) T E n F' n (0)V n (\) = 



K n -i{X, A -1 ) A'(0), j/A = A(0)^0, 
e n _ia n _iA'(0), ifX = A(0) = 0. 



Proof. Under the hypothesis, E n (t), F n (t), as well as the eigenvectors X^(X(t)) and 
X^(X(t)), are differentiable at t = (notice that E n (t) has to be constant in a 
neighborhood of t = since {aj}™ =1 C C\T). Thus, taking derivatives at t = in 
the identity 

xumfEnm^x^m) = \{t)xum) T E n {t)xi{\{t)), 

and using the fact that F n X n (X) = XX n (\) and F^ E n X n *(X) — XE n X, w (X), we 
get that 

X W (X) T E n F^(0)X n (X) = X'{0)X n *(X) T E n X n (X). 

From this equation and the proportionality between the vectors A"„(A), X„*(A) and 
V n (X), V^A), we find the desired result just noticing that 

I e„_ia„_i, if A = 0. 

□ 

Remark 4-1- From the expression (1.9) for the kernel, and taking into account that 
A is a zero of ip n , for A ^ we get that 

(4-1) ^ n _i(A,A= T ) = e n A 1 ->UA)^(A). 

Notice that y„ and c/3* can never have a common non vanishing zero since, other- 
wise, the recurrence relations (1.4) and (1.5) would imply that and ip* have this 
common zero too for all j < n, which is impossible. Thus, K n -i(X, A -1 ) = iff the 
zero A of <f> n is multiple. 

On the other hand, if A = 0, then <f> n has a zero at the origin and, so, it must 
be a n = 0. In this situation, from the recurrence relation (1.1) we have that 
4>n(z) = z<ft n ^i(z) and, hence, a n -\ = iff the zero A = of 4> n is multiple. 

As a first application of previous proposition, we will discuss the variation of 
the zeros of a n-th orthogonal polynomial under the perturbation of the last Schur 
parameter a„, which means the study of the zeros of the extensions of a finite 
segment of orthogonal polynomials. This is equivalent to analyze the eigenvalues 
of F(t) = -F n (ai, d2, . . . , a n _i, t) as a function of the complex variable t. 

Theorem 4.3. Let {<f>j}™Zo be the finite segment of orthogonal polynomials associ- 
ated to {aj}™~l C C\T. Then, the zeros of ^(z) — z<f> n -\{z) +t(f>^ l _ 1 (z) are simple 
for t e C\5, where S has only a finite number of points in any compact subset of 
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C. For t in any simply connected domain ofC\S, the zeros of (fi^ can be expressed 
as holomorphic functions oft. If X(t) is one of such functions, then, 



X'(t) = 



-e„_iA(*) J 



(^-i(AW)) 2 



if B -i(A(i) ) A(t)- 1 ) 



-e„_i- 



a n -i 



*/A(t) = 0. 



Proof. Since F(t) := T n {a\,a2, ■ ■ ■ .a„_i,i) is an analytic function of t in C, the 
number of distinct eigenvalues of F n (t) and their algebraic multiplicities are con- 
stant up to, at most, a set S with only a finite number of points in any compact 
subset of C. Moreover, for t in any simply connected domain of C\S, these eigen- 
values can be expressed as holomorphic functions A(t) [41]. 

So, just applying Proposition 4.1, and simplifying the result using Proposition 
2.3 and the relation between standard orthogonal L-polynomials and orthogonal 
polynomials, we get that 

f K n _ 1 (\(t),Mf)- r )\'(t) = -e^Ati) 1 -™^*-!^))) 2 , if A(t) ^ 0, 
I a n -i\'(t) = -1, if A(t) = 0. 

From Remark 4.1 we see that the factors of A'(t) in above equations vanish iff the 
zero A(i) of 0* is multiple. But, the right hand side can not be null (for X(t) ^ it 
can not be ip* n _ l {X{t)) = since, from (1.4), it would imply ip n _\(X{t)) — 0). So, 
we conclude that if A(io) is multiple, then A(t) is not diffcrcntiablc at t = to- From 
Theorem 3.1, the multiplicity of X(t) as a zero of </>^ coincides with its algebraic 
multiplicity as an eigenvalue of F n {t). Therefore, if F n (to) has less than n distinct 
eigenvalues, then some of them can not be expressed as diffcrcntiable functions at 
t = to and, so, t € S. On the contrary, all the points t where F n (t) has n distinct 
eigenvalues are outside S [41]. So, we find that S coincides exactly with the set of 
points where some zeros are multiple. 

Finally, notice that all the quantities with index less than n do not depend on 
t since aj is constant for j < n. So, just using similar proofs to those given in 
Proposition 2.5, Theorem 3.1 and Proposition 4.1, we see that all these results 
remain true even for t € T. □ 

Above theorem has the following immediate consequence. 

Theorem 4.4. Given a finite segment of orthogonal polynomials, the set of its 
extensions with multiple zeros is at most denumerable. Moreover, in the positive 
definite case, this set is at most finite. 

In a second application of Proposition 4.1, we will study the effect of a rotation 
of the Schur parameters on the zeros of orthogonal polynomials. First of all we will 
discuss the rotation of only one Schur parameter. This is equivalent to analyze the 
eigenvalues of F n (t) = T n {a\,a2, • • • , a/c-i, e lt ak, a,k+i, ■ • • , a n ) as a function of the 
real variable t £ [0, 2w]. 

Theorem 4.5. Let (</?^) n >o be the orthonormal polynomials with positive leading 
coefficients associated to the Schur parameters (a n (t)) n £N, where a n (t) = a n e C\T 
for n ^ k and dk(t) = e lt ak, € C\T. Then, the number of distinct zeros of ' <p„ 
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and their multiplicities are constant for t G [0, 27r]\5 7 where S is at most finite. 
The zeros of ip l n can be expressed as differentiable functions of t in [0, 27r]\5. If 
A(t) is one of such functions for n> k, then 

Kl_ 1 (X(t),M^)X\t)=-i\(t) 1 - k l [ e k . 1 a^ 

whenever A(i) ^ 0. 

Proof. First of all, notice that F n (t) = T n (a\, 02, • ■ • , fflfe-ij e lt ak, asfc+i, ■ • • , a n ) is 
not an analytic function of t considered as a complex variable. So, if we want to 
apply similar arguments to those given in the proof of previous theorem, we have 
to change the starting point. Let us consider the five-diagonal matrix T given in 
Definition 3.1 as a function of a„,a„, p n , p ni n e N, taking them as independent 
variables. That is, denoting a := (a n )„ e N an d given the arbitrary sequences in C 
b = (b n )neN, c = (c„)„ eN , d = (rf„)„ eN , let us define 



Q(a,b,c,d) := 



/ -01 Ci 
-d\a 2 —b 1 a 2 
d\d 2 bid 2 






-c 2 a 3 
-b 2 a 3 
-d 3 a 4 
d 3 di 



c 2 c 3 

b 2 c 3 
-6304 — c 4 a 5 
b 3 di -6405 



C4C5 
64C5 



For w £ C*, let Q{a{w) ) b{w) 1 c,d), where a„(w) = a„, 6„(w) = 6„ for n ^ k and 
flfc(w) = wak, bk(w) — w~ 1 bk- Then, from Q n (a(w), b(w), c, d) we can recover F n (t) 
just choosing b = a := (a TO ) ngw, c = p := (p„)„ e w, d = p := (p„)„ eN and to = e Jt , 
that is, = Q n {a{e lt ),a{e lt ), p, p). Since Gn{a(w),b(w),c 1 d) is an analytic 

function of w in C*, analogously to previous theorem we conclude that the number 
of distinct eigenvalues and their multiplicities are constant for w e C*\S", where 
S' has only a finite number of points in any compact subset of C. Moreover, these 
eigenvalues can be expressed as analytic functions of w in any simply connected 
domain of C*\S'. Therefore, the zeros of ip^, which are the eigenvalues of F n (t), arc 
constant in number and multiplicity and can be expressed as differentiable functions 
of the real variable t for t G [0, 2tt)\S, where S = S' fl [0, 2tt] must be finite. 

The rest of the theorem follows straightforward from the application of Propo- 
sition 4.1, and the simplification of the result so obtained by using Proposition 
2.3 and the relation between standard orthogonal L-polynomials and orthogonal 
polynomials. □ 

Finally, we will discuss the variation of the zeros of the orthogonal polynomials 
under the simultaneous rotation of all the Schur parameters. That is, we will study 
the eigenvalues of F n (t) = T(e tt a\,e tt a 2 , . . . , e lt a n ) as a function of the real variable 
t e [0,2tt]. 

Theorem 4.6. Let (ip n ) n —o be the orthonormal polynomials with positive leading 
coefficients associated to the Schur parameters (a n (t)) n ^, where a n (t) — e zt a n , 
a n G C\T, for all n. Then, the number of distinct zeros of ip l n and their multiplicities 
are constant fort € [0, 27r]\S' 7 where S is at most finite. Moreover, fort e [0, 27r]\5, 
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the non null zeros of p l n are simple and can be expressed as differentiable functions 
oft. If X(t) is one of such functions, then 

A'(t) = iX(t) 



K n ^(X(t),X(t)-^) 

Proof. Again, F n (t) = T{e lt a,\, e lt a 2 , ■ ■ ■ , e lt a n ) is not an analytic function of t 
considered as a complex variable. So, following the notations of previous theorem, 
we consider now, for each w £ C*, the infinite matrix Q(a(w), b(w), c, d), where 
a n {w) — wa n , b n (w) — w~ 1 b n for all n. Since F n (t) = Q n (a(e lt ) 7 a(e lt ), p, p) and 
Q n {a{w), b(w), c, d) is an analytic function of w in C*, analogously to the proof of 
previous theorem we conclude that the eigenvalues of F n (t) are constant in number 
and multiplicity and can be expressed as differentiable functions of the real variable 
t for t £ [0, 27r]\/S l , where S is at most finite. 

Let X(t) be one of such functions. Just applying Proposition 4.1 and using the 
relation between standard orthogonal L-polynomials and orthogonal polynomials, 
together with the fact that X(t) is a zero of ip^, we get that, if X(t) ^ 0, 

K£_!(A(t),A(tpj X'(t) = -i ^ai(t) + £ J i(*)j - 

where 

I^Xitr^e^a^it^fiXit))^ 

From (1.4) and (1.7) we find that 

Ij(t) - X{t)^e jPj+W ) +1 {X{t))^*{X{t)) Xitf-iej-iPj^Xit^UiXit)), 
which, taking again into account that X(t) is a zero of ip^, implies that 



K^iXit), A(i)" 1 ) X'(t) = t(pip\(X(t)) - ai(t)) = iX(t). 

From above result we see using Remark 4.1 that, if X(t n ) ^ is a multiple zero, 
then A(t) can not be differentiable at t = to- Therefore, any non vanishing zero 
X(t) must be simple for t £ [0, 2ir]\S. □ 

Remark 4.2. Concerning theorems 4.5 and 4.6, we have to remark that, since F n (t) 
is differentiable with respect to the real variable t, its eigenvalues can be expressed 
as differentiable functions A(t) in any interval of [0, 27r] where F n (t) is diagonable 
[41]. Therefore, in any interval of [0, 2tt] where the zeros of <p* n are simple, they can 
be expressed as differentiable functions. 

If r(t) = \X{t)\ and 9(t) is a differentiable determination for the phase of X(t), 
then, for X(t) ^ 0, 



X(t) r(t) 



'(*)■ 



Hence, Theorem 4.6 gives the following meaning for the n — 1-th kernel K n _i(z, y) 
associated to a sequence of orthogonal polynomials on T: if A is a non null simple 
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zero of the n-th polynomial, the real part of 1/K n _i(\, A -1 ) measures the speed of 
the rotation of A under rotation of the Schur parameters, while its imaginary part 
determines the rapidity in the radial approach of A to the origin. 

Example. As an example, we will show the consequences of this last result in the 
case of the Geronimus polynomials, defined by a constant sequence of Schur pa- 
rameters a n = a, < \a\ < 1, n > 1. In this case, it is known that the measure of 
orthogonality is supported in the arc A Q = {e l6, |o; < 8 < 2ir — a}, cos a = 1 — 2|a| 2 , 
a G [0, 7r], plus a possible mass point at zq = (1 — a)/ (1 — a) that appears iff 
3?a > |a| 2 [16, 18]. If we write z = e^, then cos tp - cos a = 2(9?a - |a| 2 ) 2 /|l - a| 2 
so, the point zq is always outside A a , except in the case 9?a = |a| 2 , for which z is 
an extremum z± = e ±la of the arc A Q . 

The orthonormal polynomials and their reversed are given by [16, 18] 

<Pn(z\ a) = 4r( M «+i ~ i 1 ~ a ) u n), 
p n 

<p* n (z; a) = y L (u n+1 - (1 - a)zu n ), 

where p = yjl — |a| 2 , u n = (iv™ — w 2 )/(wi — w 2 ) and w\,w 2 are the solutions of 
the quadratic equation w 2 — (z + l)w + p 2 z = 0, that is, 



, 2 = i{z + l±v / (^ + l) 2 - Vz}. 



Notice that w\ + w 2 = z + 1, w\w 2 = p 2 z, (wi — w 2 ) 2 = (z — z + )(z — z_) and the 
derivatives with respect to z are given by 

w[ 1 iu,-l ... 
— = , 1^3- 

Wi Z Wi — Wj 

The zeros of ip n (z; a) are the solutions of u n+ \ = (1 — a)u n , which implies 

(4 2) ^ = W2 - ^ - Q ) 

{ ' ' w% wi - (1 - a) ' 

If z is a zero of <p n (z; a), (4.1) yields 

z l — n 

(4.8) K n _i(z, z~ x ; a) = -^((1 -a)- (1 - a)z)W(u n , u n+1 ). 



being 



W(u n ,u n+1 ) =det[ U ? Un+1 



u 



n+l 



the Wronskian determinant of u n ,u n+ i considered as functions of z. From the 
definition of u n we get 

W(u n ,u n+1 )= ( 1 )2 E (-l) <+ ^«,< +1 ) 



(u>i — w 2 ) 2 w 



i w ' + _ n ( Wl - w 2 ) ( — - — J - K + w 2 ) 
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Using the properties of the functions w\, W2, together with the equation (4.2) for 
the zeros, we find finally that 

K ^ Z ~ ;a) ^ (z-z + )(z-z_) • 

Let us consider the case 3?a < |a| 2 , where the support of the measure is just the 
arc A a . Since z = \z\e te is a zero of tp n (z; a), it has to lie on the convex hull of the 
support of the measure, that is, in the subset of the open unit disk D determined 
by cos6> < cosa. Thus, \z — 1| > 1 — cosa = 2|a| 2 , \z — z \ > cosip — cosa = 
2(|a| 2 -3f?a)7|l-a| 2 and \(z - z+){z - z_)\ < 2(1 + cosa) = 4p 2 . Hence, 

(4.3) | g „-, fe -;„)| > { W - W -«■) _.}. 

Above inequality implies that, for n big enough, K n -\{z, a) can not vanish 
when z is a zero of <p„(z;a), and, so, ip n (z;a) can not have multiple zeros. More 
precisely, if JJa < |a| 2 , then 

(4.4) n > — ^|Jp — ^ + 1^ => the zeros of ip n (z; a) are simple. 

We can apply now Theorem 4.6 to the perturbation a(t) = ae lt , being a such 
that 3?a < |a| 2 . Since |a| 2 — 3?a = \a — 1/2| 2 — 1/4, we are dealing with Schur 
parameters in the region of D outside the closed disk Do = {z£ C||a — 1/2 1 < 1/2}. 
Without loss of generality we can suppose a <G (0, 1). Given < to < t\ < 2n 
such that a(to),a(ti) G ©\© , for t e [*o,ti] it must be a(t) G D\D and c(t) := 
|a(t)| 2 - Sfta(t) > c := rnin{ c(t ), c(ti)}. Therefore, from (4.4) we find that 

n > — ^ ( - — a + 1 ) =>■ the zeros of < / 9„(z; ae 2 *) are simple Vi € [to, til- 
2a 2 V c o / 

Thus, from Remark 4.2 we see that, under above condition for n, the zeros of 
(p n (z]ae zt ) can be expressed as differentiable functions z(t) in the interval [io,*i]- 
From Theorem 4.6 and inequality (4.3) we have that 



z'(t) 



z(t) 



< C n (a, to, t\) :— 



2p z 1 + a 



Co (2na 2 — l)co — (1 + a) ' 
and, thus, if z{t) = \z{t)\e l9 ^\ 9(t) differentiable in [*o,ti], we get the bounds 
||*(ti)| - |^o)||,|^i) -0(t )\ < C n (a, t , ti)|*i — * |- 

This means that, for n big enough, each zero of <p n (z; ae tto ) has a zero of ip n (z; ae ltl 
at a radial and angular distance less than C n (a, to, ii)|ti — t a \. Notice that 

, p 2 (l + a) 1 ^ / 1 
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When 3la(to),3la(ti) < we can choose a coefficient C n (a) independent of to,t\ 
because, then, cq > a 2 and, hence, 



i j. j. \ ^ s~i { \ 2 P 1 + a r l + a + a z 

C n {a,toM) < C n {a) := — — ■ — ^, for n > 



a 2 2na 4 - (1 + a + a 2 ) ' 2a 4 

Notice that C n (a) is analytic for a = 1, with C(l) = 0. This means that the best 
bounds for the location of zeros appear when a is close to 1. On the contrary, the 
behavior of C„(a) for a = is singular. This suggests a more chaotic behavior 
for the zeros under perturbations of the Schur parameters as far as the last ones 
approach to the origin, and a slower and more regular variation when they are close 
to the unit circle. 

These are just some examples of the utility of the five-diagonal matrix represen- 
tation given for orthogonal polynomials on the unit circle. A similar discussion to 
the one given here is possible for para-orthogonal polynomials, but this together 
with some applications related to the orthogonality measure will be developed in 
a separated paper [9]. There it will be discussed too an operator theoretic ap- 
proach to the study of the orthogonality measure based on the five-diagonal matrix 
representation obtained for the multiplication operator. 
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